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Abstract 

Core-collapse supernovae are connected with formation of neutron stars. Part of the gravitation 
energy is transformed into the energy of the explosion, observed in SN II, SN Ib,c type supernovae. 
The mechanism of transformation is not simple, because the overwhelming majority of the energy 
is going into weakly interacting neutrino. The attempts to use this energy for the explosion 
were not successful during about 40 years of investigation. We consider the explosion mechanism 
in which the source of energy is the rotation, and magnetic field serves for the transformation 
of the rotation energy into the energy of explosion. 2-D MHD simulations of this mechanism 
were performed. After the collapse the core consists of a rapidly rotating proto-neutron star 
with a differentially rotating envelope. The toroidal part of the magnetic energy generated by the 
differential rotation grows as quadratic function with time at the initial stage of the evolution of the 
magnetic field. The linear growth of the toroidal magnetic field is terminated by the development 
of magnetohydrodynamic instability, when the twisted toroidal component strongly exceeds the 
poloidal field, leading to a drastic acceleration in the growth of magnetic energy. At the moment 
when the magnetic pressure becomes comparable to the gas pressure at the periphery of the proto- 
neutron star the MHD compression wave appears and goes through the envelope of the collapsed 
core. It transforms into the fast MHD shock and produces a supernova explosion. Our simulations 
give the energy of the explosion 0.6 ■ 10^^ ergs. The amount of the mass ejected by the explosion 
is ~ O.14M0. The implicit numerical method, based on the Lagrangian triangular grid of variable 
structure, was used for the simulations. 
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1 Introduction 

Content. 

1. Presupernovae. 

2. Magnetorotational mechanism of explosion: 1-D calculations. 

3. 2-D MHD: Numerical method. 

4. Core collapse and formation of rapidly rotating neutron star. 

5. Magnetorotational supernova explosion. 

6. Development of magnetorotational instability. 

7. Jet formation in magnetorotational SN explosion. 

8. Mirror symmetry breaking: Rapidly moving pulsars. One-side jets. 

Supernova is one of the most powerful explosions in the Universe, producing the energy (radiation 
and kinetic) about 10^^ erg. The explosion happens at the end of the evolution of massive stars, with 
initial mass more than ^^8 Solar masses. Hertzsprung-Russell diagram of the zero age main sequence, 
and several evolutional tracks are represented in left side of Fig^ from I^Oj . Tracks in HR diagram of 
a representative selection of stars from the main sequence till the end of the evolution are given in the 
right side of Fig^ from . 

The explosion mechanisms are: 
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1. Thermonuclear explosion of C-0 degenerate core (SN la) 

2. Core collapse and formation of a neutron star, with a gravitational energy 
release ^ 6 • 10^'^ erg, carried away by neutrino (SN II, SN Ib,c). 

The transformation of the neutrino energy into kinetic one is very problematic, and was not obtained 
firmly. Magnetorotational explosion (MRE), in the core collapse SN, is based on the transformation 
of the rotational energy of the neutron star into the explosion energy by means of the magnetic field. 

Most of supernova explosions and ejections are not spherically symmetrical. A lot of stars are 
rotating and have magnetic fields. Often we can see one-side ejections. One of the main difficulties 
for core collapse supernova explosions is, how to transform any kind of energy of the star into the 
explosion energy. The magnetorotational mechanism suggested in ^ , transforms the rotational energy 
of the star into the explosion energy with a high efficiency of about 10%. It is acting in presence of 
the differential rotation when the rotational energy can be transformed to the explosion energy by the 
magnetic field action. 

2 1-D calculations 

l-D calculations of the magnetorotational explosion (MRE) have been performed first in [21 • The 
important parameter of the problem is the ratio of the magnetic and gravitational energies at the 
initial moment, when differential rotation starts to twist and to amplify the toroidal component, 
a = ^'"""^ . In ^2 rather large values of a = 0.1 and 0.01 had been considered, see Figs l2l3l 

1-D calculations of MRE for lower magnetic fields, aX a — 10^^, 10^^, 10^^, had been done in [S], 
see Fig01 Angular velocity distribution at different time moments, is represented in Fig[Slform jl2|. 
Variations of different types of kinetic energy during MRE are given in Fig|Slfrom [T^. 

The main results of 1-D calculations are the following. Magneto-rotational explosion (MRE) has an 
efficiency about 10% of rotational energy. For the neutron star with mass « I.2M0 the ejected mass 
is about O.IM0, the explosion energy is about 10^^ erg. Ejected mass and explosion energy depend 
weekly on the parameter a, explosion time strongly depends on a, Explosion time texpi ~ 1/%/^ 
|S] . Small a is difficult for numerical calculations with EXPLICIT numerical schemes because of the 
Courant restriction on the time step for a stiff set of equations: small a determines a stiffness. In 2-D 
numerical an IMPLICIT numerical schemes should be used. 

3 The difference scheme for 2-D calculations 

First 2-D calculations, showing jets in the collapse of rotating magnetized star had been performed in 
|18| . Very large magnetic field was taken, so twisting of magnetic field lines, expected in a realistic 
situation, was not important for the explosion, see Figs l7l8l 

The difference scheme used in the present calculations is described in [7], [5]. The scheme is based 
on the method of basic operators - grid analogs of the main differential operators: 

GRAD (scalar) (differential) GRAD (scalar) (grid analog) 

DlV(vector) (differential) DlV(vector) (grid analog) 

CURL (vector) (differential) CURL (vector) (grid analog) 

GRAD (vector) (differential) GRAD (vector) (grid analog) 

DlV(tensor) (differential) DlV(tensor) (grid analog) 
The scheme is implicit. It was developed, and its stability and convergence was investigated by 
the group of N.V.Ardeljan in Moscow State University. The scheme is fully conservative, providing 
conservation of the mass, momentum and total energy, and giving correct calculation of the transitions 
between different types of energies. Matrices for the transformation are symmetrical. The difference 
scheme is Lagrangian, with triangular grid, and possibility of grid reconstruction, see Fig|3 The 
angular momentum conserves automatically in this scheme. 

The astrophysical application of this scheme was started by calculations of a collapse of a non- 
magnetized rotating protostellar cloud, important in the problem of star formation Jj . Rigidly rotating 
uniform gas sphere was considered for initial conditions. The cloud, as a rigidly rotating uniform gas 
sphere, was taken with the following initial parameters: 

p = 1.492 • 10"^ Vc"^^ p = 1.548 • 10"^•'c^?/n/cm^ r = 3.81 • lO^^cm, (1) 



2 




Figure 1: Left. Hertzsprung-Russell diagram of the ZAMS and several evolutionary tracks. The masses of the 
models are indicated at the starts of the tracks, in solar units. The solid portions of the tracks indicate where 
evolution is on a relatively slow nuclear time-scale, the dotted parts show evolution on a thermal time-scale, 
and the dashed parts show an intermediate time-scale. The different symbols indicate the positions of binary 
components with well-determined masses, radii and luminosities. The position of the Sun is indicated by a 
solar symbol 0. The crosses show the position on the giant branches of the 1-, 2- and 4-Mq models where 
the binding energy of the envelope becomes positive, from )2UI . 

Right. Coarse evolutionary tracks for stars with Mi = 1, 5, 25 Mq. Heavy portions represent principal burning 
phases in the core. For M; < 2.3 Mq, CHF (core helium flash) occurs after which quiescent ^He burning begins. 
After *He has been exhausted in the core, the star arrives on the AGB (asymptotic giant branch). When the 
core with no helium in it reaches the mass ~ 0.53 Mq, TF (thermal flashes in '*He shell) start. An AGB star 
loses mass, and this process terminates with a rapid ejection of the residual hydrogen envelope in the form of 
FN (planetary nebula). CO core with M/ « 0.6 M© transforms into white dwarf. More massive AGB and 
post- AGB stars with Mi < 9 Mq evolve in a similar way, Mf increases with increasing Mi and equals 1.08 M© 
at Mi — 8.8 Mq. The symbol O indicates the onset of planetary nebula luminescence, when Te/ of the star 
reaches 3 • 10* K and the gas ionization in FN begins, from |16|. 
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Figure 2: The configuration of the magnetic field line in the subsequent moments of time for a — 0.01. The 
dashed region is the central incompressible core of the neutron star; the pointed region is a contracting part of 
the envelope. In the right figure (c) the scale is diminished to a half in comparison with (a) and (b), from |12|. 




Figure 3: The same picture as in Fig|H but for a — 0.1. The initial picture coincides with Fig|5| In the figure 
(b) the scale is diminished four times, from |12| . 




Figure 4: left A schematic picture of the initial state for the unit length of the cylinder, from |12|. 

right Shape of the field lines in the region near the core at the time t = 7/-y/a, for a = 0.01 (dashed line), and 

a = 10~* (solid line), from 
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Figure 5: The dependence of the angular velocity uj on the mass coordinate s at later time for two values of 
the parameter a. The division of the envelope in two parts is seen: the uniformly rotating central part, and 
the ejected outer part, in which the rotation became slower. The change of the direction of rotation is seen on 
the right figure, from |12|. 




Figure 6: The dependence of the kinetic energy on time for two values of a: the rotational energy of the core 
E^, envelope E^, the radial energy of the envelope Er, from 
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Figure 7: Isodensity contours in units of 10^ g/cvo? (a), and magnetic-flux contours parallel to Z-axis in units 
of 10^^ Gs/cm-2 (b), at 0.72 sec, from 




Figure 8: Isodensity contours in units of 10^ g/cvo? (a), and magnetic-flux contours parallel to Z-axis in units 
of 10^2 Gs/cm-2 (b), at 2.67 sec, from [TH|. 
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Figure 9: Elementary reconstructions: (Left up) BD connection is introduced instead of AC connection. The 
total number of knots and cells in the grid is not changed. (Right up) Adding a knot at the middle of the 
connection: the knot E is added to the existing knots ABCD in the middle of BD connection, 2 connections 
AE and EC appear, and the total number of cells is increased by 2 cells. (Down) Removing a knot: the knot 
E is removed from the grid, and the total number of cells is decreased by 2 cells, from Q 
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uj = 2.008 • IQ-^'^rad/sec, M = 1.73Mq = 3.457 • lO^^g, 7 = 5/3, u'' ^ ^ 0. 

On the outer boundary the pressure was taken to be equal to a small constant {p = 0.87 • 10~^^ 
dyn/cm^). At the outer boundary of the cloud, the gravitational potential $ is defined by the integral 
Poisson formula, using the expression for the volume potential. The results of calculations, showing 
formation of the bounce shock, and subsequent outburst are represented in Fig^|from ilj. 

4 Presupernova Core Collapse 

Collapse and formation of a rapidly and differentially rotating neutron star had been calculated in [2] 
in non- magnetized approximation. Equations of state takes into account degeneracy of electrons and 
neutrons, relativity for the electrons, nuclear transitions and nuclear interactions. Temperature effects 
were taken into account approximately by the addition of the pressure of radiation and of an ideal gas. 
Neutrino losses and iron dissociation were taken into account in the energy equations. A cool white 
dwarf was considered at the stability limit with a mass equal to the Chandrasekhar limit. To obtain 
the collapse we increase the density at each point by 20% and switch on a uniform rotation. 

The hydrodynamical equations with gravity for modelling the nonstationary processes in rotating 
gaseous bodies are: 

^ = v, ^+pV.v = 0, = -V (P) - pV$, (2) 

+ PV • V + pF{p, T) = 0. A$ = A-nGp. 

at 

Here the equation of state and the function of neutrino losses have been represented by following 
functions 

P = P(p,T) = Po(p)+p5RT+^, e = eo{p) + l^T+ — +eFe{p.T), (3) 

6 I p 

where epe{p, T) is the iron dissociation energy, F(p,T) is the function of neutrino losses, among which 
only the URCA process losses are important. Other losses, such as pair annihilation, photo production 
of neutrino, plasma neutrino were also included in the calculations. 



P^'^ = hp^'^l{l + c,p^/% p 



< 



\ P(|'=' = a . l0&'^to-8-4i9)='= < p < Pk. ^' 

p 

-5Rr+ +eFe{p,T), £o(p) = / — r^— dp, £Fe{p,T) = 

2 p J p-' Amp \TiFe-ToFe 







^6 



1.3 • 10^si(T)T" 1 , 

/(P,T) = ' erg.g-i-s-\ ^ = T.10-^ (6) 

1 + (7.1 • 10-5p/T )2/5 

r _ 1, T_<7, 

ffi(T) = <^ 664.31 + 51.024(T - 20), 7 < T < 20, (7) 
[ 664.31, r>20. 

The URCA losses f{p,T) have been described by formula |I2|, obtained by approximation of tables 
from J7]. Other losses (unimportant) are included in Qtot- 

F{p,T) = {f{p,T) + Qtot)e'^. (8) 

Neutrino absorption is important in deeper layers of the neutron star. It is taken into account 
implicitly, by multiplying the transparent neutrino flux by the multiplier e~To , where the effective 
neutrino optical depth Ti, is calculated using local density gradients 12 Initial state was chosen 
as a spherically symmetric star, with a mass 20% larger than the limiting mass of the corresponding 
white dwarf M = l,0042Mo + 20%, and rotating uniformly, with an angular velocity 2.519 (1/sec). 
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Figure 10: Time evolution of the shape of the rotating collapsing gas cloud, from 
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Figure 11: The triangle grid at the beginning, t = 0, from 



The temperature distribution was taken in the form T = bp^l^ . The initial triangular grid covering a 
quarter of the circle, where equations are solved, is represented in Fig llll 

Results of calculations are represented in Figs^]-^^ The collapse is stopped due to formation of 
a stable neutron core, rotating almost uniformly. The shock is reflected from the core, leaving behind 
a hot plasma, where vortexes are generated. It leads to a strong mixing, what is important for the 
interpretation of observations of the light curve of SN 1987A. The light curve, corresponding to the 
radioactive decay of '^^Ni — > ^^Co ^^Fe was observed in this object, what could be explained by 
enrichment of the outbursting matter due to mixing with deeper layers. 

The main result of these calculations is a self-consistent model of the rapidly rotating neutron star 
with a differentially rotating envelope. The adjusting parameters are chosen to reproduce the calcula- 
tions with a refined account of the neutrino transport jl3| . In correspondence with these calculations 
we have obtained that the bounce shock wave and neutrino deposition do not produce SN explosion. 
Distribution of the angular velocity is represented in Fig^] The period of rotation of the almost 
uniformly rotating core of the young neutron star is about 0.001 sec. 



5 2-D magnet or ot at ional supernova 

Calculations of magnetorotational core-collapse supernova have been performed in J5J . Magnetohydro- 
dynamic (MHD) equations with self-gravitation, and infinite conductivity have been solved using the 
same numerical scheme as described above. The set of equations is the following: 



H 



dx dp ^ . dv ,/ H-H 

d^=^' d^ + ^^-^ = «' ^dF--^™H'^+^ 

de 



V ■ (H H) 

47r 



p-(-)=H.Vv, A<I> = 4^Gp, p-+FV-v + pi;^(p,T) = 0, P = P(p, T), e = £(p, T), (9) 



At 
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Figure 12: Density distribution in the central region at maximum compression state t — 0.1425 sec, from 




Figure 13: Velocity field, showing mixing of matter behind the shock at i = 0.1828, from 
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Figure 14: Angular velocity distribution in the equatorial plane of the new born neutron star at t — 0.261 
after the beginning of the collapse. 

where ^ = ^ + v • V is the total time derivative, x = (r, z), v = {vr, v^,Vz) is the velocity vector, p 
is the density, P is the pressure, H — {Hr, H^, Hz) is the magnetic field vector, $ is the gravitational 
potential, e is the internal energy, G is gravitational constant, H (g) H is the tensor of rank 2, and 
F{p, T) is the rate of neutrino losses. Equation of state and the function of neutrino losses have been 
the same as in the previous section. Additional condition here is divH = 0. The problem has an axial 
symmetry ( ^ =0), and the symmetry to the equatorial plane (z=0). Boundary conditions are the 
following: 

P = p = T = Bcf, = at the outer boundary; Vr = j,- = i?r = at r = 0; (10) 

dB 

— — — a.t r — 0; Vz — 0, — (dipole-like) or = (quadrupole-like) at z = 0. 

Initial toroidal current (see Figll5|l was taken at the initial moment (time started now from the 
stationary rotating neutron star) producing , according to Bio-Savara law 

V 

Initial magnetic field of quadrupole-like symmetry is obtained by opposite directions of the current 
in both hemispheres (Fig ll6() . Magnetic field is amplified due to twisting by the differential rotation, 
and subsequent development of the magnetorotational instability. The field distribution for initial 
quadrupole-like magnetic field with a = 10~^, at the moment of the maximal energy of the toroidal 
magnetic field is represented in Fig ll7l Two dark areas: near the equatorial plane, and around the 
axis at 45°, show the regions with local maxima of the toroidal magnetic field H^. The maximal 
value of Hff, = 2.5 • 10"'^^ Gs was obtained in the calculations. The magnetic field at the surface of the 
neutron star after the explosion is ff = 4 • 10^^ Gs. Temperature and velocity fields development 
with time is given in Fig^l Time dependencies during the explosion of different types of the energy: 
rotational energy, gravitational energy, internal energy, kinetic poloidal energy, are given in Figs ll9l 
The evolution of the magnetic poloidal energy and magnetic toroidal energies are given in the right 
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R 

Figure 15: Initial, at t = 0, current distribution in the upper hemisphere. Dipole-like field is obtained when 
currents in both hemispheres have the same sign, and quadrupole-like field appears when currents have opposite 
signs. 




Figure 16: Initial quadrupole-like magnetic field, at f = 0. 
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Figure 17: Toroidal magnetic field distribution at the moment of its maximal energy. 

side of Fig l2()l Distribution of the specific angular momentum rv^ at different time moments is given 
in FigEOKleft), and FigEI] 

Almost all gravitational energy, transforming into heat during the collapse, is carried away by 
weakly interacting neutrino. Dependence of the neutrino luminosity, and integral neutrino losses on 
time are represented in Fig|221 A mass "particle" is considered ejected if its kinetic energy is greater 
than its potential energy. Time dependence of the ejected mass in (Mq), and the ejected energy 
during the magnetorotational explosion with a quadrupole-like field are represented in Fig l23l The 
total energy ejected in the kinetic form is equal to 0.6 • 10^^ erg, and the total ejected mass is equal to 
0.14M„. 



6 Magnetorotational instability 

Magnetorotational instability (MRI) leads to exponential growth of magnetic fields. Different types 
of MRI have been studied in JS] , [22] , [1] , ■ MRI starts to develop when the ratio of the toroidal 
to the poloidal magnetic energies is becoming large. In 1-D calculations MRI is absent because of a 
restricted degree of freedom. Therefore the time of the MR explosion is increasing with a as iexpi ^ , 

a = . Due to the development of MRI the time of the MR explosion depends on a much weaker. 

The ratio of two magnetic energies is changing with time almost with the same speed for all a, so MRI 
starts almost at the same time. The MR explosion happens when the magnetic energy is becoming 
comparable to the internal energy, at least in some parts of the star. While the starting magnetic energy 
linearly depends on a, and MRI leads to exponential growth of the magnetic energy, the total time of 
MRE in 2-D is growing logarithmically with decreasing of a, texpi ~ — log a. These dependencies 
are seen clearly from 1-D and 2-D calculations with different a giving the following explosion times 
^Gxpi (in arbitrary units): 

a 0.01, icxpi = 10, a = 10"^^ tcxpi = 10*^ in 1 - D (12) 
a = 10~^ iexpi --6, a = 10-^^, iexpi -12 in 2 - D 
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Figure 18: Time evolution of the velocity field (left column) and temperature (right column) for the time 
moments t = 0.07s, 0.20s, 0.30s, from g]. 
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Figure 19: Time dependence of the gravitational, internal, rotational, and kinetic poloidal energies of the star 
during magnetorotational explosion with a quadrupole-like field, from _4,. 
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Figure 20: Right Time dependencies of the magnetic (toroidal and poloidal separately) energies of the star 
during magnetorotational explosion with a quadrupole-like field. Left The specific angular momentum v^pV 
distribution for the time moment t = 0.07s, from |1] 
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Figure 21: Time evolution of the specific angular momentum for the time moments t — 0.20s, 0.30s. The 
darker parts of the plots correspond to the higher specific angular momentum, from 
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Figure 23: Time dependence of the ejected mass in (Mq), and ejected energy during the magnetorotational 
explosion with a quadrupole-like field. The particle is considered as ejected if its kinetic energy is greater than 
its potential energy and its velocity vector is directed from the center, from 




Figure 24: Duration of the explosion time on a in presence of MRI, from |19| 
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6.1 Toy model of the MRI development 



Toy model of the MRI development shows an exponential growth of the magnetic fields at initial stages 
of MRI development in 2-D. Initially MRI leads to formation of multiple poloidal differentially rotating 
vortexes. Angular velocity of vortexes is growing (linearly) with a growth of H^. 

The right-hand side is constant at the initial stage of the process. When the toroidal field reaches its 
critical value if*, the MRI instability starts to develop. As follows from our calculations, the critical 
value corresponds to the ratio between total over the star toroidal and poloidal magnetic energies 
~ 1 — 3. A local value of this critical ratio is much larger, reaching at the beginning of MRI ~ 10'' in 
the regions of the maximal growth of the toroidal magnetic field, what in time roughly corresponds to 
about 100 rotations of the neutron star core. These regions serve as germs for subsequent development 
of MRI in other parts of the star. The appearance of MRI is characterized by formation of multiple 
poloidal differentially rotating vortexes, which twist the initial poloidal field leading to its amplification 
according to 

where I is the coordinate, directed along the vortex radius, w,y is the angular velocity of the poloidal 
vortex. Qualitatively the poloidal field amplification due to the vortexes induced by MRI is shown 
in the Fig. [23 The enhanced poloidal field immediately starts to take part in the toroidal field 
amplification according to (|13|l . With further growing of H^p the poloidal vortex speed increases. Our 
calculations give the values of = 0.0132 s"! at \H^ \ = 2.46 • lO^^G corresponding to t = 0.041 s and 
ujv = 0.052 s""' at \Hip\ = 4.25 • lO^^G corresponding to t — 0.052 s^^ for the same Lagrangian particle. 
In general we may approximate the value in brackets 114|l by linear function on the value {H^ — i/*) 
a,s 

^l^=a{H^-H;). (15) 
^{H^-H;) =AHroaiH^~H;), (16) 
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7 Jet formation in MRE 

Jet formation in MRE happens when the initial magnetic field is of a dipole-like structure. 2-D 
calculations with the initial dipole-like magnetic field gave almost the same values of the energy of 
explosion ~ 0.5 • 10^^ erg, and ejected mass « O.14M0 (see Figs l^ . but the outburst was shghtly 
collimated along the rotational axis |19j . 

7.1 Violation of mirror symmetry of magnetic field 

There is a violation of the mirror symmetry of magnetic field when two components dipole -f quadrupole 
are present from the beginning |23| . Another possibility of violation of the mirror symmetry of magnetic 
field was considered in '111 . If we start from axially symmetric configuration with a mixture of poloidal 
and toroidal fields with different parity, than any exchange between these two components will lead 
to such violation. For example, if we start from the dipole-like poloidal field (negative parity), and 
symmetric toroidal field with a positive parity, than twisting of the dipole magnetic field by the 
differential rotation leads to creation of additional toroidal field with a negative parity. Therefore the 
total resulting toroidal field will become asymmetric relative to the equatorial plane, see Fig l27l 
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Figure 25: Qualitative picture of the development of the MRI in 2-D. 
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Figure 26: (Left) Time dependence of the ejected energy during the magnetorotational explosion with a 
dipole-like field, (right) Time dependence of the ejected mass during the magnetorotational explosion with a 
dipole-like field, from |19) . 
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Figure 27: Loss of a mirror symmetry of the magnetic field in the differentially rotating star, from 

The magnetorotational explosion with a mirror asymmetric field will be always asymmetrical. A 
kick velocity due to such asymmetry of the magnetic field could be up to ~ 300 km/sec The 
interaction of the neutrino with asymmetric magnetic field could increase the kick velocity up to ~ 1000 
km/sec, giving a possibility to explain the origin of velocities of the most rapidly moving pulsars |1()| . 

8 Conclusions 



1. In the magnetorotational explosion (MRE) the efficiency of transformation of rotational energy 
into the energy of explosion is 10%. This is enough for producing core collapse SN from rapidly 
rotating magnetized neutron star. 

2. Development of magneto-rotational instability strongly accelerates MRE, at lower values of the 
initial magnetic fields. 

3. The new born neutron star has inside a large (about 10^^ Gauss) chaotic magnetic field. 

4. Jet formation is possible for dipole-like initial topology of the field, what may have a possible 
relation to cosmic gamma-ray bursts; equatorial ejection happens at prevailing of the quadrupole-like 
component. 

5. Braking of the equatorial symmetry happens in differentially rotating star with toroidal and 
poloidal components of an opposite symmetry; it may explain formation of the rapidly moving pulsars 
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and one-side jet formation. 
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